Stark effect of shallow impurities in Si 
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We have theoretically studied the effect of an electric field on the energy levels of shallow donors 
and acceptors in silicon. An analysis of the electric field dependence of the lowest energy states in 
donors and acceptors is presented, taking the bandstructure into account. A description as hydrogen- 
like impurities was used for accurate computation of energy levels and lifetimes up to large (several 
MV/m) electric fields. All results are discussed in connection with atomic scale electronics and solid 
state quantum computation. 
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PACS numbers: 71.70.Ej,71.55.Cn,03.65.Fd,03.67.Lx 



I. INTRODUCTION 

The field of atomic scale electronics (ASE) aims at con- 
trolling charge and spin in semiconductors at the level of 
individual dopant atoms. Such an ability is very attrac- 
tive, both for physics and for the development of (quan- 
tum) devices. From a fundamental point of view, dopant 
atoms are interesting, because they can be considered as 
the solid state analogue of atoms in free space. Several 
well-known effects from atomic physics (e.g. the Stark ef- 
fect and Zeeman effect) have been studied in great detail 
in large ensembles of dopant atoms 1 . The prospect of 
experimentally realizing atomic scale electronics has re- 
newed the interest in dopant atoms. Measurement and 
control of individual dopant atoms allows for the observa- 
tion of quantum coherent time evolution and interactions 
of the dopant's wave functions, which is essential for the 
operation of a quantum computer. 

Manipulation of a single particle's wave functions can 
be realized by using a local magnetic or electric field. 
Such a field can be used either to perform the desired 
manipulation itself, or to provide a local perturbation 
allowing for addressing a single impurity by a global ra- 
diation field. A local electric field could be realized by 
putting a small gate close to a dopant atom, which is 
in principle accomplishable with current technology. An 
ultimate application of gate-manipulation is found in the 
solid state quantum computer as proposed by Kane 2 ^. 

To get more insight in the physics of atomic scale elec- 
tronic devices, it is essential to try to predict their poten- 
tial behavior. A first step is the description of isolated 
dopant atoms in a (homogeneous) electric field. Much 
more difficult is accurate modelling of a the time evolu- 
tion of a dopant atom wave function in an inhomogeneous 
field and the description of the interaction of two or more 
dopants in a field. 

Dopant atoms binding one electron or hole can be de- 
scribed as a hydrogen atom, where the vacuum values of 
the dielectric constant and the electron mass are replaced 
by the appropriate values for the semiconductor. This 
'scaled hydrogen model' (SHM) provides a reasonable de- 
scription of the dopant atom's energy levels. Therefore, 
it is useful to look at existing studies of the Stark ef- 
fect in the hydrogen atom. Calculation of the shift and 



splitting of the hydrogen energy levels up to very large 
electric field have been carried out by several different 
method a 4 i 5 i 6 . Within the SHM, these results can be di- 
rectly translated to dopant atoms in a uniform electric 
field. However, we found that almost no actual results of 
such calculations in the (field)range of interest for ASE 
have been published. 

The SHM also offers a manageable way to describe 
a dopant atom in an inhomogeneous electric field. Re- 
cently, several calculations using this framework have 
been published£*&il i n the context of quantum comput- 
ing. However, the SHM fails in the explanation of ef- 
fects where it is essential that the bandstructure of the 
semiconductor is taken into account (as an example, see 
Ref.E3). 

Many measurements of the energy levels of dopant 
atoms in semiconductors (large ensembles) are known, 
but only a few concerning the effect of a uniform elec- 
tric field have been reported, presumably because such 
measurements are much more difficult than e.g. mea- 
surements in a magnetic field or under stress. Among 
them are spectroscopic measurements of the boron en- 
ergy levels in silicon subject to electric fields up to 
0.15 MV/mii. Electron-spin-resonance experiments 12 
demonstrated that the electric field couples linearly to 
the acceptor ground state. The magnitude of the effec- 
tive electric dipole moment for linear Stark coupling has 
been estimated as 0.26 D for boron acceptors in silicon 
(1 D = 3.3 x 10~ 30 Cm). Photo-ionization measure- 
ments have shown a very large electric field effect on the 
phosphorus ground state in Si 1 ^, but this was measured 
in highly doped samples where the interaction between 
dopants dominates the Stark effect of individual energy 
levels. Finally, quadratic level shifts have been observed 
in deep selenium double donors in Si, located in the space 
charge region of a diodeA^. 

In this paper, we will theoretically investigate the ef- 
fect of a uniform electric field on isolated shallow impu- 
rities in silicon. Primary interest for ASE will be in the 
ground state and possibly the first few excited states. 
These states are the only ones that are well separated 
from neighboring levels and at low temperatures only 
the ground state is occupied. Therefore, we focus on 
the lowest energy states of impurities in silicon. First, 
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we derive the shift, splitting and wave functions of the 
lowest donor levels in silicon in a small uniform electric 
field, taking full account of the multiple valley conduc- 
tion band structure (SeclBj. We briefly outline a similar 
calculation for acceptors in silicon fSec. The results 
are useful for applications where a local gate is used to 
bring a single dopant atom into resonance with a global 
radiation field (nuclear magnetic resonance, electron spin 
resonance). Moreover, they can be used to outline the 
limitations of the SHM. Second, in Sec. |W] we present 
accurate numerical calculation of the Stark effect in sil- 
icon within SHM, from zero field up to fields that are 
relevant for atomic scale electronics and quantum com- 
puting (several MV/m; see for instance Ref. 0). Finally, 
we conclude by discussing possible extensions and alter- 
natives of our methods which are useful to address issues 
in ASE (Sec. El . 



II. DONORS 




FIG. 1: (a) Schematic representation of the conduction band 
valleys of silicon as constant energy surfaces in fc-space. The 
six valleys are labelled by numbers, e.g. 4 represents the [010] 
valley, (b) Definition of the coordinate system with respect to 
the Si-crystal unit cell. We have x || [100], y || [010], z || [001], 
v || [110], and w || [111]. The orientation of the figure in part 
(a) and (b) is the same. 



A. The donor ground state 

Group theory is a powerful tool to derive various prop- 
erties of dopant wave functions in a semiconductor. In 
order to provide the necessary background and to fix the 
notation, we will briefly review some relevant properties 
of donor levels in silicon (see e.g. Ref. 0). Degeneracy 
due to spin is not lifted by an electric field in donors. For 
simplicity, we will therefore not count those degeneracies 
in this section. 

The conduction band of silicon has six equivalent min- 
ima located on the [100] and equivalent axes. These min- 
ima are commonly called 'valleys' and we label them by 
the numbers 1 to 6 as shown in Fig. ^a). The band 
structure in the vicinity of valley 1 , located in /c-space at 
ki = (fco,0, 0), can be approximated as 
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where mu = 0.98m and m± = 0.19m are the electron ef- 
fective masses and m is the free electron mass. Further- 
more, ko = 0.85— 16 >, where a is the size of the silicon 
unit cell. Similar expressions hold for the remaining five 
valleys. 

From effective mass theory (EMT) it follows that the 
ground state wave function of the Hamiltonian of an elec- 
tron bound to a donor can be written aa— 



$(r) 



6 

£ 



(1) 



where the are numerical coefficients and the F^fr) are 
envelope wave functions, which are slowly varying on the 
length scale of a. Fi(r) = F^r) satisfy the hydrogen- like 



Schrodinger equation 
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F(r) = EF(r) (2) 



and similar equations hold for the remaining Fp. The 
ipp (r) are Bloch-wave functions at the minimum of val- 
ley fi and can be written as e 4k ^' r it M (r), where ^(r) has 
the periodicity of the silicon crystal lattice. Because for 
all /j the eigenvalues resulting from Eq. J5J) are the same, 
Eq. shows that the degeneracy of each of these eigen- 
values is multiplied by six for the total wave functions 
^(r). In particular, the ground state solution of Eq. @ 
gives rise to a six-fold degenerate donor ground state. 

The symmetry group of the conduction band minima 
(and thus of the Bloch functions <p(r)) is Coq V in EMT, 
which reduces to Civ in the silicon crystal. The envelope 
wave functions F(r) belong to -Dooh- Their products be- 
long to the cross-section of both groups, which is C^v 
For the ls-like (m = 0) ground state function of Eq. |2 
F)i(r), such a product transforms according to the Ti 
representation of the valley symmetry group C^v Be- 
cause the donor is located at a substitutional site of the 
tctrahcdral silicon lattice, the total wave function has T^- 
symmetry. Using Frobenius' theorem^, it can be shown 
that the Ti representation of Ci v induces the Ti +T 3 + T 5 
representation^ of T^. This means that linear combina- 
tions of the F)j(r) can be found that have the correct 
transformations properties under T ( i- Using the notation 
a = (ai, . . . , ag) (as in Eq. the reduction to the Td 
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representations is carried out by 

a g = ^=(1,1,1,1,1,1) T, 

o»r = ^(-1,-1,-1,-1,2,2) 1 

a s = 1(1, 1,-1,-1,0,0) J 3 

(3) 

£** = ^(1,-1,0,0,0,0) ^ 
a y = ^(0,0,1,-1,0,0) I r 5 

a z = ^(0,0,0,0,1,-1) J 

Each of the vectors a defines a wave function \I> through 
Eq. ^ Here, the basis functions of the two- and three 
dimensional representations have been chosen such that 
& r and \I/ S transform under Td as 2>z 2 — r 2 and V3(x 2 — 
y 2 ), respectively. Similarly, $> x , ^> y and ty z have been 
chosen such that they transform under Td as x, y and z, 
respectively. 

The potential term in the EMT-Schrodinger equation 
@ is a good approximation only for r > a, where a is 
the lattice constant of silicon. For small r, the charge 
of the nucleus is not screened by other electrons and it 
will attract electrons much stronger than described by 
the potential in Eq. J2J. Because the symmetry of the 
potential is not affected, the states are still described by 
the representations of Td, but they are no longer degen- 
erate. The Ti state ^ g is the only one of the six ground 
state wave functions that has non-zero electron density 
at the nucleus (r = 0). Therefore, it has a larger binding 
energy than predicted by EMT and for most donors in 
silicon the ls(ri) state is the true ground state. This 
effect is generally called 'chemical splitting' (because the 
size of the effect depends on the donor in question) or 
'valley-orbit splitting'. The remaining states (especially 
the non-s states) are quite well described by EMT, be- 
cause the electron density at the nucleus is negligible. As 
an example, in case of phosphorus in silicon, the ls(Ti) 
state (the ground state) has been measured to be located 
45.29 meV below the conduction band minimumi, while 
the EMT-prediction is 31.27 meVii 

B. Symmetry of the donor ground state in an 
electric field 

After this brief review of established knowledge of sili- 
con donors, we return to the main subject of this paper. 
From purely symmetry based considerations, we can find 
how the Hilbert subspace spanned by the original six val- 
ley wave functions is decomposed by the application of an 
electric field in a certain direction. The impurities con- 
sidered in this paper occupy substitutional sites in the 
silicon lattice and their wave functions transform accord- 
ing to representations of site symmetry group Td. The 
symmetry group of a uniform electric field £ is Coc V . 
When £ is applied in an arbitrary direction in the sili- 
con crystal, the symmetry group Td of the Hamiltonian 
reduces to the trivial group C\. Only when the direction 



TABLE I: Reduction of the site symmetry of an impurity in 
a uniform electric field in various directions and the resulting 
reduction of the irreducible representations—. 



Direction 


(100) 


(111) 


(110) 


Group 




Cs v 


c s 


Ti (T d ) 




Ti 


ri 


r 2 (T d ) 


r 3 


r. 


r 2 


r 3 (T d ) 


Ti + r 3 


r 3 


r a +r 2 


r 4 (T d ) 


r 2 + r 3 + r 4 


r 2 + r 3 


ri + 2r 2 


r 5 (T d ) 


ri + r 2 + r 4 


ri +r 3 


2Ti + r 2 


r 6 (f d ) 


r 5 


r 4 


r 3 + 4 


r 7 (T d ) 


r 5 


r 4 


r 3 +4 


r 8 (Td) 


2T B 


r 4 + r 5+ 6 


2r3 +4 



of the field is along one of the main crystallographic di- 
rections of the crystal, the result is Ci v for £ \\ (100), 
C 3v for £ || (111), and C s for £ || (110). The reduction 
of symmetry can induce a splitting in the original energy 
levels as shown in Table [I] As expected, the electric field 
does not remove degeneracy due to time reversal symme- 
try and therefore all resulting levels are at least two-fold 
degenerate. 

To make the connection to the valley wave functions 
T) 1 (r)0 AI (r), we will now discuss the symmetry of the Is 
levels in an electric field from another point of view. We 
start by looking at the individual valley wave functions 
and subsequently derive which linear combinations form 
appropriate donor wave functions (using the method of 
Ref. |21|) . When a donor impurity in silicon is situated 
in an electric field along the positive z-direction, the val- 
leys 5 and 6 keep their Ci v symmetry, while the field 
reduces the symmetry group of the other four valleys to 
C\ . These four valleys are mixed by the elements of the 
site symmetry group Ci v and are therefore grouped to- 
gether in the third column of Table [H] 

In case of a Is state, the valley wave functions belong 
to the Ti representation of Civ (for valley 5 and 6) or C\ 
(for 1, 2, 3 and 4). This is found by reducing the even 
m = representation of Dooh to Ci v and C\, respec- 
tively. By using Frobenius' theorem, it can be deduced 
that these generate for the impurity wave function the 
representations Ti and ri + Ti + T 3 + T4 of Ci v , respec- 
tively. This is also shown in Table ITT1 together with the 
(set of) induced wave function(s) spanning the subspace 
of that representation. In a similar way, we obtained 
results for the electric field in the other main crystallo- 
graphic directions. They are also shown in the table. 

Due to the valley-orbit splitting (which has been ig- 
nored so far) the three irreducible components of the 
donor ground state are already energetically separated 
at zero field. Therefore, the basis vectors have to be 
chosen in such a way that they agree with the zero-field 
energy splitting of the Ti, T 3 and T 5 levels of T^. The 
result for various directions of the electric field is shown 
in Table Hill 
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TABLE II: Considering the symmetry of the valley wave func- 
tions in an electric field, the symmetry of the total wave func- 
tion they induce can be obtained. The results for the Is 
level, without considering valley-orbit splitting, are shown in 
this table. The direction of £ in the first column is denoted 
by the vectors defined in Fig. ^a). The fifth column lists 
the representations of the appropriate site symmetry group, 
given in the second column. The basis vectors are given in 
the notation of Eq. J3J . 

Dir. Site Valley Valley E(site) Basis 
£ sym. sym. 



z C21, 1, 2, 


3, 4 d 




(1,1,1,1,0,0) 






r 2 


(1,-1,1,-1,0,0) 






r 3 


(1,1,-1,-1,0,0) 






r 4 


(1,-1,-1,1,0,0) 


5 




Ei 


(0,0,0,0,1,0) 


6 




Ei 


(0,0,0,0,0,1) 


w C 3v 1, 3, 


5 C a 


Ei 


(1,0,1,0,1,0) 






E 3 


(lo 2 , 0,w,0, 1,0) 








(w,0, Lo 2 ., 0,1,0) 


2, 3, 


6 C s 


Ei 


(0,1,0,1,0,1) 






E 3 


(0, LO 2 , 0,LJ, 0,1) 








(0,o;,0,^ 2 ,0,l) 


v C s 1, 3 


Ci 


Ei 


(1,0,1,0,0,0) 






r 2 


(1,0,-1,0,0,0) 


2,4 


Ci 


Ei 


(0,1,0,1,0,0) 






r 2 


(0,1,0,-1,0,0) 


5 


C s 


Ei 


(0,0,0,0,1,0) 


6 


Cs 


Ei 


(0,0,0,0,0,1) 



C. Shift and splitting in an electric field 

Now, we will derive the shift and splitting of the low- 
est donor levels in an electric field from a perturbation 
calculation. Results for other levels can be derived using 
the same method, although (because the level spacing is 
smaller for higher levels) the range of fields where the 
perturbation calculation is valid is much smaller. 

Although the six- fold degeneracy of the ls-levels is 
lifted by the valley orbit interaction, the complete man- 
ifold is relatively well-separated from the higher levels 
(the separation of the highest Is ^3) level to closest ex- 
ited level (2po) is roughly twice as large as the separation 
between the ls(ri) and ls(r 3 ) levels). Therefore, we 
consider the ls-manifold as a whole in a single perturba- 
tion calculation, taking only the coupling among the Is 
levels themselves into account. 

The electric field couples to the (induced) dipolc mo- 
ment D = er of the impurity state and gives rise to an 
additional term in its Hamiltonian — £ ■ D, reflecting the 
energy associated with the dipole in the field. By mak- 



ing use of the Wigner-Eckart orthogonality theorem from 
group theory2£, it is possible to find the vanishing ma- 
trix elements as well as the dependencies between the 
non- vanishing matrix elements, as they follow from the 
symmetry of the system. The Is sub- matrix [H] of the 
total Stark Hamiltonian H = Ho + 5 • D is given by 



PVoZy 











E 3 






E?, 



-pib£x P3b^/^£x 
-p3$£y p35V~3£y 
2p 35 £ z 



Pl5^x 



P\^£y Pv£z \ 



~P35£x 
P35V3£x 



-PZ^y 



%P35£z 




E 5 

Pb£z 

P5^y 



Pb£z 

E 5 
p$£x 



P5£y 
Pf>£x 
E 5 



The elements of this matrix are given by \H\ij = 
((pi\H.\tpj ), where the wave functions ifi are taken from 
the basis (* g , <3> r , * z )) as defined before. The 

energies E\ , E3 and E§ are the eigenvalues of the unper- 
turbed Hamiltonian 7io, that is the zero-field energies 
of the ls(ri), ls(F3) and ls(r5) level, respectively. For 
phosphorous in silicon, the values are E\ = —45.59 meV, 
E 3 = -32.58 meV, and E 5 = -33.89 meV with respect 
to the conduction band edgei. The parameters P15, P35 
and j»5 describe the coupling between the ls-levels. As 
can be seen, these are the only three independent parame- 
ters describing the coupling between the levels. They can 
be expressed in terms of integrals over products of wave 
functions, e.g. P15 = e(^l' g |a;|^' 2 ;) and P5 



Perturbation theory is invoked by calculating the 
eigenvalues and eigenvectors of this 6x6 matrix up to 
second order in £. This yields the Is energy levels and 
wave functions as a function of electric field for £ along 
the three main crystallographic directions. The energy 
levels are presented in the last column of Table ITTT1 From 
Table iHfl it can be seen that the ls(Ei) ground state ex- 
periences an isotropic quadratic shift downwards^, while 
for the other levels the behavior depends on the direction 
of the electric field. In Figure|21the results for £ \\ (100) 
are plotted schematically. 

The corresponding eigenvectors were also obtained 
from this calculation. In the limit £ — > they coincide 
with the vectors given in Table ITTT1 allowing to label each 
eigenvalue with the correct representation. These results 
are directly applicable in the prediction of allowed optical 
transitions between the various levels. 

We discuss the behavior of the three Is states in in 
some more detail. The normalized eigenfunctions in an 
electric field parallel to z (again up to second order in 
£) corresponding to the eigenvalues already given in Ta- 
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TABLE III: Reduction of the Is donor energy levels in an electric field. The basis vectors belonging to these states are given 
(in the notation of Eq. ©) in the limit £ — > {uj = e 27rz ^ 3 ). The eigenvalues (up to second order in £ ) are the result of the 
perturbation calculation described in the text. 



Field direction £ = £ 7^ Basis vector(s) 



Eigenvalue 



ri(Td) ri(c 2 „) (1,1,1,1,1,1)7^6 



Ex 



r 3 (T d ) ri(<7 ao ) 
r 8 (C a „) 

r 5 (T d ) ri(c a .) 
r 2 (c 2 „) 



l,l,l,l,-2,-2)/Vl2 
1, 1,-1,1, 0,0) /2 

0. 0.0.0.1,-1)/V2 

1, -1,1,-1, 0,0) /V2 
1, -1,-1,1, 0,0) /V5 



£3 + 
£3 



|2p 35 r 



£0+ |p 5 |£ 

E 5 - \p 5 \£ 



)S 2 



Ti(T d ) ri(c 3 „) 
r 3 (T d ) r 3 (c 3 „) 



r 5 (T d ) ri(c 3 „) 
r 3 (c 3 „) 



1,1,1,1,1,1)/V6 

a; 2 , a; 2 , w, w, 1, l)/\/6 
a;,tj,w 2 ,tj 2 , 1, 1)/V6 



-£ 2 



Si 



IP 1 2|f>35 1 C-2 
£ 3 + B^TEjt 



1,-1,1,-1,1,-1)/V6 £ 5 ±|^3|p 5 |£: + ( 
aj 2 ,-w 2 ,w,-a;, 1,-1)/V6 £ 5 =F |V3|psl£ 
w, — a>, w 2 , — cj 2 , 1, — 1)/V6 



P15l" _ 4|p 35 | 
E 5 -E 1 E3--E5 



)£ 2 



Ti(T d ) ri(C7.) 

r 3 (T d ) riCC) 
r 2 (c s ) 

r 5 (T d ) riCC) 
ri(c s ) 
r 2 (c fl ) 



1, 1, 1, 1, 1, 1)/V6 £1- 

1, 1, 1, 1, -2, -2)/a/12 £ 3 + 

1, 1,-1,1, 0,0) /2 E 3 + 

0. 0.0.0.1,-1)/V2 

1, -1,1,-1, 0,0) /V2 £5 
1, -1,-1,1, 0,0) A/2 £ 5 



B5--B1 



P35 



B3--B1' 

3|P35| 2 

E3--B1' 



£5 + \p B \£ - U-E 3 - El 



\p35\_ 



\P3bL 



L£15l S |C 2 

IpisI 2 \c2 

EE-Si 



)£ 2 



E3—E1 




Field 



bleCTDare 



where 



- iA/2(* x + * a ) 
= i72(* x - *„) 
$ 2 = /3£ • - 0'£ ■ 9 r - 



Pl5 



E 5 — Ei 



(l-i(|/3| a + |/3?)£ 2 )¥, 



2 P35 , = h 2 P35 



E% — En 



E3 — Ei 



(4) 



FIG. 2: Schematic plot of the Is energy levels as a function of 
the electric field £ . The values of the parameters ps, pig, and 
p 3 5 have been chosen such that the plot clearly illustrates the 
qualitative features of the Stark effect in the energy levels. 



The initial zero field wave function ^ g has the highest 
spacial symmetry possible in a tetrahedral lattice. To 
get more insight in the contribution of the six valleys as 
a function of the applied field, we can write the perturbed 
ground state wave function $ 9 in the notation of Eq. © 
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as 

(1, 1,1, 1,1,1) + (o,o, 0,0,-y, i')£+ 

(-7, -7, -7, -7, -7', -i)£ 2 , 

where 

7 = ^(I/?| 2 +/3'V2), 7 ' = i(|/?| 2 -2/3"V2), 7 " = /3V3 

and an overall factor 1 / \/6 was omitted. From these 
expressions, we see that the contribution of the valley in 
the — z direction increases linearly with the field, while 
contribution of the opposite valley decreases linearly with 
the field. This reflects the field-induced dipole moment 
of the ground state. 

The results of this calculation could be made quan- 
titative if the values of the parameters ps, P15, and P35 
were known. This can be done by evaluating the integrals 
defining these parameters and using e.g. the EMT wave 
functions from Eq. Q • However, due to the strongly os- 
cillating integrants, this is numerically a non-trivial task. 
Furthermore, the EMT wave functions have a higher 
symmetry than the lattice, and the value for p$ obtained 
in this way is always zero. An estimate for p§ can only 
be obtained using more sophisticated approximations for 
the wave functions. More importantly, the applicability 
of such results is limited, especially for the Is state, as 
the effects of valley-orbit interaction are not included in 
the EMT wave functions. 

It is important to note that the energies in Table IIIII 
and the eigenstates in Eq. are based on symmetry 
properties only and not on the explicit form of the EMT 
wave functions. Therefore, these results remain valid, 
even if valley-orbit interaction and central cell correc- 
tions are fully included. Such modifications would only 
influence the values of the parameters p$, P15, and P35. 

III. ACCEPTORS 

Acceptor wave functions can be equally well used for 
ASE as donors. Recent experiments showing that the co- 
herence time of spins of bound holes is more than 1 ms2£, 
even justify the prospective use of acceptor wave func- 
tions as qubits. We therefore also briefly outline the prop- 
erties of silicon acceptors in an electric field, taking the 
silicon valence band structure into account. The initially 
threefold degenerate valence band maximum is split by 
spin orbit interaction, which causes one of the bands to 
shift downwards by ~ 43 meVi Due to the spin-orbit 
interaction, spin is not a good quantum number anymore 
and the bands must be characterized by the total angular 
momentum, which is § for the upper two bands. Due to 
the half- valued angular momentum, the Bloch wave func- 
tion at the valence band maximum transforms according 
to one of the double valued representations of T^, namely 
r§. As a result, the total impurity wave functions trans- 
form according to representations of the same group. The 



ground state wave function, as well as the first few ex- 
cited levels belong to the r§ representation and they are 
all four- fold degenerate (including spin). 

A. Linear Stark effect 

To derive the small-field splitting of acceptors in silicon 
in an electric field, we use degenerate perturbation the- 
ory for each level individually. To that end, the Hamil- 
tonian sub-matrix (tpi\H.\(pj) of the level under consider- 
ation must be calculated and diagonalizcd, where the ipi 
form a suitable basis for the subspace of that particular 
level. 

As mentioned before, the components of the electric 
dipole operator er transform according to the rows of the 
T5 representation of T^. Because the anti-symmetrized 
direct products {r 6 x r 6 } = {r 7 x r 7 } = Ti do not 
contain T^, the first order Stark matrix vanishes for lev- 
els with Tq or Tj symmetry. Hence, such levels do not 
experience a linear Stark effect. On the other hand, 
{r 8 x r 8 } = Ti + T 3 + T 5 does contain r 5 , so that a 
linear Stark effect is possible for a Tg level^. 

The effective linear Stark Hamiltonian 34 for a Tg level 
is given by 24 

[Wig 111 = ^ P8 (S X {Jy, J Z } + £y{J Zl J X } + £ Z {J X , Jy}), 

where the parameter pg is related to the effective dipole 
moment of such a state. The Ji (i = x, y, z) are matri- 
ces of the components of the angular momentum opera- 
tor with respect to some convenient basis and {A, B} = 
i (AB + BA) is the anti-commutator. The eigenvalues of 
this matrix are given by 

e 8 ±\ P8 \£, 

where both eigenvalues occur twice. This is a symmetric 
splitting of the level, which is independent of the direc- 
tion of £. Note that p 8 vanishes within EMT, similar to 
P5 before. Estimates of p 8 obtained in literature range 
from 1(T 2 DS4 to 0.26 JM. 



B. Quadratic Stark effect 

Because {r 6 x r 6 } = {T 7 x r 7 } = r l5 the quadratic 
effective Stark-Hamiltonian for a Tq and T 7 level is simply 
given by 

^cff,quad — a i£ ^ i 

where / is the identity matrix and the (i — 6, 7) are 
phenomenological parameters, that can be expressed in 
terms of integrals over wave functions. It follows that the 
Tq and T 7 levels experience an isotropic quadratic shift 

E t +a z £ 2 , 
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where E{ is the unperturbed energy of a Ti level. The 
two-fold degeneracy due to time reversal symmetry is ob- 
viously not removed by the electric field. 

The quadratic part of the effective Hamiltonian for a 
Tg level, such as the ground state, is given by 2 ^ 

[H]r d = a s £ 2 I + b 8 [J%£l + J 2 £ 2 + J 2 £ 2 - \ J 2 ] 

""^"^73^8 [{tZc, Jy\£ x £y "h \_Jyi J z\£y£z ~t~ Jx\£z£x\ j 

where as , 6s and Cg are again phenomenological parame- 
ters. The total Hamiltonian has two distinct eigenvalues 

a 8 £ 2 ± [p 2 £ 2 + b 2 £ 4 

+ (c 2 - 3b 2 )(£ 2 £ 2 + £ 2 £ 2 + £ 2 £ 2 ) 

+ 6p$c$£ x £ y £ z ] 1/2 , 

(5) 

each of which is still doubly degenerate (due to time re- 
versal symmetry) 35 . For £ | (100) this expression re- 
duces to (up to second order in £) 

E s ± \ Ps \£ + a s £ 2 . 

For £ || (111) we find 

E 8 ±\p 8 \£ + (a s ±^V3c 8 )£ 2 

and for £ || (110) we have 

E & ± \p s \£ + a 8 £ 2 . 

The results for £ || (100) and for £ || (110) are the same 
in this approximation, but different in third order. 

Obviously, the wave functions of donors and acceptors 
are very different and this is reflected in their respective 
electric field behavior. The donor ground state under- 
goes an isotropic quadratic shift. The acceptor ground 
state has an isotropic linear splitting, superposed on an 
anisotropic quadratic shift. 

In the spectroscopic measurements of boron acceptors 
in siliconil (mentioned in the introduction), the observed 
Fg-levels indeed show a quadratic shift. However, the 
expected level-splitting was not observed, most likely due 
to limited resolution. 



IV. LARGE ELECTRIC FIELDS IN SHM 

In this section, we will calculate energy levels of an 
impurity in a semiconductor as a function of electric field 
in the range from zero to ~ 5 MV/m. This is done within 
the scaled hydrogen model, where the band structure of 
the semiconductor is accounted for by a single effective 
mass and the dielectric constant only. 

For this calculation it is convenient to express all quan- 
tities in so-called effective atomic units. For instance, 



TABLE IV: Atomic units for some relevant physical quantities 
in vacuum and silicon. For silicon the values e s = 11.4 and 
m* = 0.26 (appropriate for electrons) were taken. 



Quantity 


Unit 


Value in vacuum 


Value in Si 


Energy 


2Ry 


27.2 eV 


54 meV 


Length 


do 


0.053 nm 


2.3 nm 


Electric field 


2Ry/eao 


510 GV/m 


24 MV/m 


Time 


ft/2Ry 


2.4- 10~ 17 s 


1.2 • 10" 14 s 



energies are expressed in units of twice the effective ion- 
ization energy and length in units of the effective Bohr- 
radius. Conversion of units of relevant quantities for both 
vacuum and silicon are given in Table llVl 

In the past, several algorithms have been described 
in literature to calculate electric field dependence of the 
energy levels of the hydrogen atom. However, very little 
results in the range of interest for ASE (fields up to ~ 
0.1 a.u. 3 ) have been published. Therefore, we found it 
important to fill this gap by fully presenting the results 
of our calculation. For this purpose, we used the slightly 
adapted version of a variational algorithm that not only 
yields the energy levels, but also their lifetimes^. 

For completeness, we will very briefly outline the main 
features of this method. The hydrogen Schrodinger equa- 
tion (including the electric field) in parabolic coordinates 
can be separated, which allows for high numerical accu- 
racy without too much computational effort. In order to 
be able to find the energy positions of the resonances as 
well as their lifetimes, the complex scaling method was 
applied 25 . Then, for each coordinate the Hamiltonian 
(including electric field) is expanded with respect to a 
truncated basis of unperturbed wave functions. This can 
be done analytically. Finally, the energy levels and life- 
times are obtained by tracking (separately for each level) 
the eigenvalues of this matrix from zero field in small 
steps to larger fields. 

By using the method described above, we calculated 
the energies of all states with n = 1,2,3 for < £ < 
0.2 a.u. The results for the energy levels are depicted 
in Fig. |3| The levels are labelled by parabolic quantum 
numbers— (ni,n,2,m), which are more suitable for hy- 
drogen in an electric field than the more common spher- 
ical quantum numbers (n, l,m). The magnetic quantum 
number m has the same meaning in both representa- 
tions. The main quantum number n is related to the 
parabolic quantum numbers by n = rii + n% + |to| + 1. 
The electric field lifts all degeneracies except for spin and 
(ni, ri2, ±m). So (including spin) there are both two- fold 
degenerate levels (m — 0) and fourfold degenerate levels 
(m ^ 0). 

Figure |21 shows that the ground state (n = 1) exhibits 
a small second order shift downwards. The n ~ 2-level 
splits into three levels. Two of them are (for small £) lin- 
early shifting upwards and downwards. The middle one 
has no first order shift, consistent with the well known 
results from perturbation theory2&. Finally, the ninefold 
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£ (a.u.) 

FIG. 3: Evolution of the lowest lying energy levels (n = 1,2,3) 
of a hydrogen- like system versus electric field £. For conver- 
sion of a.u. to conventional units, see Table HVI 




£ (a.u.) 

FIG. 4: Energy width and lifetime of the lowest lying energy 
levels of hydrogen-like systems (n = 1,2,3) versus electric 
field £ . For conversion of a.u. to conventional units, see Ta- 
ble EV| 

0i 1 



degenerate n = 3-level can be seen to split into six levels. 
As expected, the effect of the electric field on higher lev- 
els is larger, due to their larger spatial extent. At large 
values of the field, several levels cross each othe»2& and 
some of them show non-monotonous behavior. 

The few results of calculations that can be found in 
literature (obtained by different methods) and overlap 
with our results are in very good agreement, both for the 
ground stated and for the first excited state (m = 

The method we used for our calculations can not only 
be extended to very large fields, but it also has the ad- 
vantage of yielding the width of the energy levels. The 
increasing energy width of the hydrogen-like levels in an 
electric field is the results of the ability of the field to 
ionize the atom. The finite probability for the carrier 
to tunnel out of the nucleus' potential well leads to a 
finite lifetime^ of the level. In Fig. 0] the evolution of 
the width of several hydrogen energy levels is depicted. 
Obviously, the width of all levels is zero at zero field, 
which is equivalent to an infinitely long lifetime. For any 
nonzero £, the lifetimes have a finite value, that decreases 
monotonously with the field. The stronger the binding 
energy of a level at zero field, the faster the lifetime de- 
creases when the field increases. 

In Figure [3 the results of Figure |3] and 0] are com- 
bined into one 'intensity map', where the levels are dis- 
played as normalized Lorentzian line shapes, the width 
of which is taken from Fig. 0] The figure shows clearly 
that for the realistic electric field £ = 0.04 a.u. (about 
1 MV/m; see Table llV)l the energy width of all levels 
except the ground state is already larger than or compa- 
rable to their binding energy. The ground state lifetime is 
only 10 ns at that field. We also note that for our purpose 
it is not very useful to extend the calculation to higher 
fields, as already at £ = 0.2 a.u. all levels are very much 
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-0.2 
S -0.3 
-0.4 
-0.5 



10" 3 10" 2 10" 1 

£ (a.u.) 

FIG. 5: Map of the energy levels from Figure [3] converted to 
Lorentzians using the data of Figure |1] For conversion of a.u. 
to conventional units, see Table llVl 



broadened and strongly overlapping. Although in case 
of hydrogen atoms in vacuum such large field (0.2 a.u. 
~ 100 GV/m) are only realized in astronomy, in semi- 
conductors they can be easily achieved under laboratory 
conditions (0.2 a.u. ~ 5 MV/m). 

Though the SHM oversimplifies the bandstructure, it 
is in our opinion particularly useful to estimate lifetimes. 
Fig.^Jshows that the lifetimes are primarily a function of 
the zero-field binding energies. Assuming this is still true 
when the silicon bandstructure is included, interpolation 
of the results can be expected to provide a good first order 
approximation of the level's true lifetime. For example, 
the n = 1 value in Fig. ^underestimates the phosphorous 
donor ground state lifetime, because it is stronger bound 
than assumed in EMT. 
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When the electric field is generated by a small local 
gate, this gate is usually separated from the semicon- 
ductor by a potential barrier that is sufficiently high to 
prevent tunneling. If the distance of the dopant atom to 
the barrier is not too small, ionization of the dopant atom 
can still occur in large fields (and the lifetimes discussed 
before still apply). However, the charge carrier will not 
be 'lost', but transferred to the potential well created by 
the biased gate&. 

V. DISCUSSION AND CONCLUSION 

In the preceding sections, we have used two distinct 
approaches to study the behavior of impurity wave func- 
tions in an electric field. The first includes details of the 
bandstructure, but is only valid for small fields and is 
somewhat qualitative. From this symmetry-based anal- 
ysis, we derived the energy level shift and splitting for 
donors and acceptors in small electric fields, as well as 
the modification of the donor wave function. Further- 
more, the symmetry classification of the resulting levels 
provides for straightforward prediction of allowed optical 
transitions. 

The second approach, the scaled hydrogen model, is 
fully quantitative and applicable up to large fields, but 
neglects most features of the silicon bandstructure. Still, 
the SHM offers a manageable and valuable way to de- 
scribe important phenomena in atomic scale electronics. 
We presented the energy width and lifetime of the impu- 
rity levels in large electric fields, calculated within this 
framework. 

It is possible to combine the two approaches and treat 
Eq. J21 in a way similar to that presented in Section llVl 
Though this is in principle straightforward, the reduced 
symmetry and lack of separability will make this ap- 
proach numerically very involved. Furthermore, it is im- 
portant to note that the direction of the electric field 



with respect to the valley axis is not the same for all val- 
leys. As an example, for £ || z the energy levels of F5 
and Fq are affected in a different way than those of the 
other four F M . If the solutions for the various valley wave 
functions are known, they can be combined into impurity 
wave functions using the data in Table ITT1 

Though potentially interesting, such an effort is not 
likely to yield a good description of the dopant's wave 
function at high electric fields, despite the tremendous 
increase of necessary computational power. The reason 
is the omission of valley orbit interaction, which not only 
affects the ground state, but also the coupling to excited 
states. Especially for large fields, the coupling influences 
the properties all energy levels. It has been shown that 
inter-valley coupling accounts for the splitting of the Is 
state for P in Si21. Inclusion of this effect appears to be 
a minimum requirement for obtaining accurate quantita- 
tive results valid at large fields. 

Recently, calculations of a silicon donor in an electric 
field in the tight binding approach have been presented^. 
This approach seems to be a useful alternative to calcula- 
tions based on effective mass theory. Given the fact that 
this method inherently includes the band structure of the 
semiconductor host, it is striking how similar the results 
are to calculations based on the SHMS. This underlines 
the power of the SHM in this type of calculations. 

In summary, we have calculated the Stark effect of im- 
purities in silicon in two different approaches. Moreover, 
we discussed the results and the computation methods 
used in the context of atomic scale electronics and quan- 
tum computation. 
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In contrast to our treatment of donors, we will use the 
technique of effective Hamiltonians to derive the matrices 
for acceptor levels. 

Note that there is a mistake in the corresponding expres- 
sion in Ref. I24I where the last term between the square 
brackets is missing. 

Levels belonging to the same representation of the spatial 
symmetry group Coov can be seen to cross each other in 
Fig-El This is however no violation of the non-crossing rule, 
since for this specific problem there exists an additional 
constant of motion that is associated with the separability 
of the Hamiltonian--. 

This lifetime is solely due to the possibility of ionization 
and is unrelated to (radiative or non-radiative) transitions 
from an excited level to a lower state. 



APPENDIX A: CHARACTER TABLES 

In this appendix, we give for completeness the charac- 
ter tables of various symmetry groups, relevant for this 
paper. Table Ivl refers to the lattice symmetry group T4. 
Depending on the direction of the electric field, it reduces 
to one of the groups Civ , C% v or C s , the character tables 
of which are given in Table ED Finally, the table of the 
continuous groups Dooft, and Coc V are given in Table IvTTl 
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TABLE V: Character table for the double group T d . 





E 


E 


8C* 3 


8C3 


3C 2 , 3<5 2 


6S4 


6S4 


6crd, 6(j£j 


ri 


1 


1 


1 


1 


1 


1 


1 


1 


r 2 


1 


1 


1 


1 


1 


-1 


-1 


-1 


1 3 


o 
A 


2 


— 1 


1 

— 1 


2 


U 


U 


U 


r 4 


3 


3 








-1 


1 


1 


-1 


r 5 


3 


3 








-1 


-1 


-1 


1 


r 6 


2 


-2 


1 


-1 





-V2 


V2 





r 7 


2 


-2 


1 


-1 





V2 


-72 





r 8 


4 


-4 


-1 


1 















TABLE VI: Character tables of the single valued irreducible 
representations of the point groups C 2t ,, Cs v , and C s . 



C2v 


E 


C2 


(Tv 


a v 


C'iv 


E 2C 3 3cr„ 


Cs 


E a 


Ti 


1 


1 


1 


1 


Li 


1 1 1 




1 1 


r 2 


1 


-1 


1 


-1 


r 2 


11-1 


r 2 


1 -1 


r 3 


1 


1 


-1 


-1 


r 3 


2-10 






r 4 


1 


-1 


-1 


1 











TABLE VII: Character table of the groups Doofe and Coou 
(upper left part). These continuous groups have four (two for 
Coov) one- dimensional representation and an infinite number 
of two-dimensional representations. Here m is a positive in- 
teger. 





E 


2C£ 


00a v 


i 


2S-£ 


00C2 


r+ 


1 


1 


1 


1 


1 


1 


rj 


1 


1 


-1 


1 


1 


-1 


■pm 
1 9 


2 


2 cos mip 





2 


2 cos mip 





r+ 


1 


1 


1 


-1 


-1 


-1 




1 


1 


-1 


-1 


-1 


1 


pm 


2 


2 cos mip 





-2 


—2 cos mip 






